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We consider thermodynamic and transport properties of a long granular array with strongly 
connected grains (inter-grain conductance g ^> 1.) We find that the system exhibits activated 
behavior of conductance and thermodynamic density of states ~ exp{— T* /T} where the gap, T* , 
is parametrically larger than the energy at which conventional perturbation theory breaks down. 
The scale T* represents energy needed to create a long single-electron charge soliton propagating 
through the array. 

PACS numbers: 73.23.-b, 73.23.Hk, 71.45. Lr, 71.30.+h 



The low-temperature conductivity of granular materi- 
als continues to attract attention, from both experimen- 
talists |H and theorists 0, El ■ From a conceptual point of 
view, an attractive feature of these systems is the possi- 
bility to separately control the effects of electron interac- 
tion and quantum interference. A particularly interesting 
situation is realized in arrays with large inter-granular 
conductance, j » 1 (in units of e 2 //i), and large grain 
size (small electron mean level spacing, S, in the grains). 
In the limit gS -C T electron transport in such systems 
becomes purely inelastic and long range quantum coher- 
ence is inhibited Q. As we show below, under these 
conditions interaction effects alone lead to an exponen- 
tial suppression of conductivity, which is fully amenable 
to analytical treatment. 

At high temperatures, the conductivity of a granu- 
lar array is Ohmic, a — g (hereinafter the length of 
the system is measured in the number of grains). At 
lower temperatures Altshuler-Aronov interaction correc- 
tions Q begin to impede the conduction behavior. For 
"inelastic" arrays this correction was found Q to be 
5a = — In E c /T, where E c is the charging energy of the 
individual grains. Comparison with the Ohmic contri- 
bution shows that this perturbative correction is small as 
long as T > E c = E c e~ 9 . At the same energy scale, E c , a 
single grain connected to external leads would cross over 
to the strong Coulomb blockade regime 0, IE ■ 

In this paper we show that the conductivity of a Id 
array of grains crosses over to a manifestly insulating 
(activated) behavior at a parametrically larger tempera- 
ture, T* ^> E c . Below the crossover, the conductivity is 
exponentially small: 

o- = 9 cxp , T < T* , (1) 

as characteristic for insulators. The size of the gap, T*, is 
model-dependent. For arrays with vanishing background 
charge, q, at each grain, we find T* ~ E c e^ g ^, while in 
the case of random background charges T* ~ E c e~ g / 2 . 
In either case T* ^> E c . In the case of q = 0, the ther- 



modynamic density of states (DOS) is suppressed along 
with the conductivity. Note that Eq. JJJ is not a result 
of phonon-mediated hopping, but is a consequence of in- 
teractions between electrons only. 

The reason why the scale T* and Eq. were over- 
looked previously is that they are not visible in standard 
perturbative expansions in l/g -C 1. In the conventional 
formulation of the theory in terms of voltage fluctua- 
tions 0, Eq. (TJ comes from including large, topolog- 
ically non-trivial fluctuations (instantons). Proliferation 
of instantons leads to insulating behavior at temperatures 
T* ^S> E c , where Gaussian fluctuations are still small. 
Notice that for a single grain instantons affect the con- 
ductance only at much lower temperatures T ~ E c |IJ. 
However, contrary to a single dot, an extended array pro- 
vides a large 'entropic volume' for the formation of in- 
stantons, which substantially increases the characteristic 
temperature. We shall return to a quantitative discussion 
of this picture below. 

It turns out, however, that the effect is more naturally 
explained using a language of charge fluctuations. It is 
known that even a highly conducting barrier retains some 
ability to pin the charge on a single grain 0, 0, Q ■ This 
mechanism is drastically enhanced in the array geome- 
try, where it bears similarity to the pinning of charge 
density waves by a periodic potential. The elementary 
mobile excitations in this system are extended solitons 
of unit charge. Their activation energy, T* , is given by 
the geometric mean of the pinning strength and inverse 
charge compressibility (grain capacitance) . Our main re- 
sult, Eq. QJ, simply reflects the thermal density of such 
charge solitons. 

To quantify this latter picture we consider a general- 
ization of a model previously employed to study quantum 
dots [Tjj. Its simplest version treats the grains coupled 
by a single conducting channel and therefore has g < 1. 
(We shall show later that the predictions derived from it 
survive generalization to the complementary case j> 1.) 
The model is formulated in terms of a charge displace- 
ment field, 0j(t), where 9j+i — 9j = Nj is the charge on 
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the j-th grain. In the absence of backscattering at the 
contacts, the action reads 

Af— l 

So=EtE[ E = - °i - if + , (2) 

3=1 m 

where the first term represents the charging energy of 
the grains, while the second originates from integrat- 
ing out the continuum of the electronic degrees of free- 
dom. Backscattering at the inter-granular junctions is 
described 0] a nonlinear term to be added to the ac- 
tion ©: S hs = ^r!2jJ drcaa(2ir9j(r)). Here, r is the 
reflection amplitude and D 3> E c the effective bandwidth 
of the model. 

A crucial observation that makes the problem solvable 
is that even for r = the quantum fluctuations of 9j (r) 
do not diverge in the limit T — > 0: 
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where = AE C sin 2 (7rfc/2M) is the excitation spectrum 
defined by Eq. (J2J), and C « 0.577 is the Euler constant. 
One can thus safely integrate out these fluctuations, to 
arrive at a sine-Gordon type action that involves only the 
classical (zero Matsubara) component of the field: 



S[6] 



El 

T 



M-l 



g ) 2 -2 7 cos(27r^)] , (4) 



3 = 1 



where 7 = \r\e c / (2tt 2 ) . In the multi-channel case the 
coupling constant generalizes to 7 ~ ]1 |r s |, where 

r s is the reflection coefficient of the sth channel. 

Equation is known as the action of the Frenkel- 
Kontorova model This model describes a harmonic 
elastic chain of "atoms" with stiffness E c , placed on top 
of a periodic "substrate" potential with the amplitude 
2~fE c . The " incommensurability parameter" q represents 
the periodicity mismatch between the chain and the sub- 
strate. For small values of q the system will find it favor- 
able to retain a commensurate state (cf. Fig. n a ))> i- e - 
the chain will stretch a little so as to still benefit from 



an optimal coupling to the substrate. Setting 9j 



0, 



one finds that the energy per atom in this configuration 
is given by F c — E c (q 2 — 2j). At \q\ > 27 this energy 
becomes positive, and the state with 0j = can not per- 
sist as the lowest-energy state (e.g. it is obviously less 
favorable than the incommensurate state with 9j jq 
and Fi w 0). Indeed, in the limit of weak periodic po- 
tential, the transition between the commensurate, (Fig. 
n a) and incommensurate (Fig. ^ b) phases occurs at 
\q\ = q* = y/2j. For the average number of electrons per 
grain, N(q) = q~d q F/ (2E C ), one thus expects: N(q) = 
for \q\ < q* (insulator) and N — ► q for \q\ > q* (metal). 



FIG. 1: Plot of the function (see text) N(q) of an atomic chain 
in proximity to a periodic substrate, a) (b)) commensurate 
(incommensurate) configuration, c) solitary excitation with 
its excitation energy (inset.) 



The relevant thermal excitations superimposed on the 
commensurate ground state are so-called incommensura- 
tions (solitons in the language of the sine-Gordon model) 
- localized defects, where the distortion 9j 'climbs' over 
a maximum of the substrate potential to relax back into 
a next minimum (cf. Fig. n°))- Minimizing the action 
subject to the boundary condition = 0(1), and 
using the condition 7 « 1, one finds that the action of 
one of these solitons is given by S s — T*(q)/T, where 



T*{q) kl = q ' 2nE c (q* 



2 7 . (5) 



As a result, at \q\ < q* the thermodynamic density of 
states scales as d q N{q) ~ exp(-T*(q)/T), i.e. T*(q) 
is the excitation gap of the system. Consequently the 
conductivity exhibits the same activation behavior, cf. 
Eq. Q). Notice, that for \q\ = q* the gap vanishes. In 
agreement with our earlier estimate, this signals a pro- 
liferation of solitary excitations and the proximity of the 
incommensurate phase. 

A more thorough discussion of the system (cf. 
Ref. 0|) shows that insulating 'plateaus' along with 
superimposed solitary excitations form not only around 
q = 0, but also around other rational values of q. How- 
ever, both the width of these plateaus and the corre- 
sponding activation energies decrease for higher rational 
fractions. Among the low lying rationals, q — 1/2 plays 
a particularly interesting role. Indeed, for a single grain, 
q = 1/2 represents charge degeneracy point, where the 
system is in a conducting state (Coulomb blockade peak). 
Unexpectedly, the array exhibits a very different behav- 
ior. Using our current language, q = ±1/2 is special in 
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that the atoms of the unperturbed chain alternatingly 
find themselves in minima/maxima of the substrate po- 
tential. Under these conditions, energy can be gained 
by building up an 'Peierls-distortion' of periodicity 2 and 
modulation amplitude 50 j ~ 7. This configuration is in- 
ert against small variations in q (insulating) where, how- 
ever, the width of the insulating plateau estimates to only 
A<7i/2 ~ 7, i.e. it is much smaller than Aq = 2q* ~ y^y. 

The above discussion was based on the arguably arti- 
ficial assumption that the background charges in every 
grain are the same. Under realistic conditions, though, 
one expects q — ► qj to fluctuate. (The same applies to 
the tunneling conductances and charging energies; we be- 
lieve, however, that these latter fluctuations are of lesser 
relevance.) Let us briefly consider the extreme limit 
where qj € [0, 1[ on the different grains are uniformly dis- 
tributed statistically independent random variables. For 
an undistorted chain, dj+i — 6j = qj, the potential terms 
27 cos 27t#j vary randomly and the energy per atom is 
zero on average. The system can then gain an energy 
SF w -E c 7 2 per grain by slightly distorting the chain, so 
that that S9j + i — 268 "j + SOj^i — j8m.(2wjqj). In anal- 
ogy with the "clean" case, the excitation energy of this 
deformed state is expected to be T* ~ ^/E C SF « E c j. 
(Note that this mechanism closely related to the collec- 
tive pinning of Abrikosov lattices in type II superconduc- 
tors 13.) 

Having discussed the charge pinning mechanism in the 
context of the few channel model, we next turn to the 
generalization to highly conducting arrays (j > 1). To 
this end we employ the so-called Ambegoakar-Eckern- 
Schon (AES) model Q. This formalism describes the 
system in terms of the quantum phase, 4>j(r), conju- 
gated to the charge 8j + i(r) — 0j(t) of the j-th grain. 
(Alternatively, one may think of <pj as the time inte- 
gral of the voltage on the grains, i(f>j — Vj.) The action 
of the model contains two terms, S = S c + St, where 
s c[4>] = S dT [4> 2 ]/{^ E c) - iqfa], is the charging energy 
of a grain kept at voltage Vj — i(j>j , and 
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describes the process of tunneling. Here Aipj = (cj>j+i — 
4>j)/2 where i<f>o and i<f>M are the voltages on the leads 
connected to the array. 

Before analyzing the array in terms of the above ac- 
tion, let us review a few general features of the AES ap- 
proach: (i) ignoring effects of quantum interference, the 
applicability of the model is restricted Q to temperatures 
T > g5; (ii) The quadratic approximation to the action, 



, provides a 
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complete description of the classical RC-resistor network 
corresponding to the array; (iii) anharmonic fluctuations 



of the phase lead to the perturbative logarithmic correc- 
tion to the dc conductivity mentioned in the introduc- 
tion; (iv) technically, the field <pj represents a mapping 
S 1 — ► S 1 from the unit circle (imaginary time augmented 
with periodic boundary conditions) into itself (4>j is a 
phase). In addition to (f> = 0, the tunneling action St[4>] 
of a single grain (which, for low temperatures T <C E c 
represents a good approximation to the total action of 
the grain) possesses a set of topologically non-trivial ex- 
tremal phase configurations known as Korshunov instan- 
tons [111HII3: 
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Here, W £ Z \ is the winding number of the mapping 
(j)^ and z = (zi, . . . , Z\w\) is a set of \W\ complex pa- 
rameters constrained by \z a \ < 1. The action associated 
with the instanton, S[cj)^] w g\W\ — 2niqW, is nearly 
z-independent (lif which identifies the z a 's as instanton 
zero modes. (Physically, argz Q determines the instance 
and 1 — \zj\ the duration of the voltage pulse, i(j)( z \) 

Turning to the array, the fact that the tunnel- 
ing action depends only on the differences of neigh- 
boring phases, A0j, implies that a 'plateau' formed 
by L instanton fields embedded into M — L zeros, 
(0, . . . , 0, 4>( z \ ■ ■ ■ , (j^ z \ 0, . . . , 0), represents an extremal 
configuration. For W = ±1 its action is given by 
S[<p] — L(tt 2 T / 1 E c ^f2wiq) + g , an expression that suggests 
an alternative interpretation of the instanton plateau: 
rather than monitoring a state of every grain, one may 
think of the plateau as a dipole of two charges located 
at the positions of the step-wise changes in the winding 
number, Wj\ — » 1 and 1 — ► 0, resp. Within this pic- 
ture, exp(— g/2) represents the fugacity of the charges, 
\L\tt T / E c , their interaction, and the g-dependent term 
describes the interaction of the dipole with a uniform 
electric field 2iriq. More formally, a summation over 
all instanton configurations followed by integration over 
massive Gaussian fluctuations and zero modes 0, 0] 
leads to the expression 
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where 7 2 = g 3 e 9 , and the interaction potentials are 
ir 2 T 

V(ja~jb) = -^e a e b \ja-jb\; $tf a ) = 2mqe a j a , (9) 

with e a = (—1)°. These equations generalize from a sin- 
gle dipole to the statistical mechanics of a Id Coulomb 
gas in a uniform external field. The fugacity of the gas, 
~/E c /T, results from multiplication of the instanton ac- 
tion by the fluctuation factor [19j. 
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To understand the properties of this system, we recall 
the standard mapp ing of a Coulomb gas onto the sine- 
Gordon model 20] . In the present context, the action of 
the latter is given by Eq. @, which completes the proof 
of equivalence of the two approaches discussed in this 
Letter upon the proper identification of 7 |2l| . Therefore 
the activation temperature, T*(q), is given by Eq. (J5J 
with 7 = g 3 / 2 e~ g / 2 . It is worthwhile to mention that a 
key element in establishing this equivalence is the factor 
E c /T ^ 1 in the fugacity of the Coulomb gas, Eq. ijHJl. 
This factor results from the large volume available to 
fluctuations in the array geometry, i.e. no such factor 
exists for single grains. 

We finally turn to the discussion of the low tempera- 
ture (T < T*) dc transport properties of the array. As 
mentioned above, in the insulating phase the fundamen- 
tal excitations of our system are solitary configurations 
carrying unit charge. Referring for detailed discussion to 
Ref. [18| . we here merely mention that in the presence 
of an external field, E, the dynamics of these objects is 
controlled by the Langevin equation 
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= gE + £{t), (10) 



and the linear conductivity exhibit activation behavior. 
The mechanism discussed here may play a important 
role in the construction of an "extended" theory encom- 
passing both strong interaction and effects of long range 
quantum interference. 

We are indebted to A. I. Larkin and J. Meyer for 
valuable discussions. This work is supported by NSF 
grants DMR97-31756, DMR01-20702, DMR02-37296, 
and EIA02-10736. 
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where is a Gaussian correlated noise with 

(£(W)>€ = 9T5(t - f ) and d%/df = 9 j+x - 20 3 + 
0j-\ is the discrete second derivative. 

In the commensurate phase < q*) the solutions of 
this equation are solitary configurations, Qj(t) = 9(j—vt), 
propagating with a constant velocity, v. Substituting 
this ansatz into Eq. (|10fl . one finds v = 7~ 1 / 2 gi?, where [11 
.y-1/2 ^ 1 ig the soliton length. As each of these ob- 
jects carries unit charge, the current density is given 
by J = env, where n = 7 1 / 2 e~ T l T is the concentra- 
tion of the thermally excited solitons (T* is given by 
Eq. (JSJ with 7 = g 3 / 2 e -9 / 2 ). The linear dc conductiv- 
ity of the array is thus given by Eq. In the case 
of a random background charge, qj, we expect a sim- 
ilar result with, however, a different activation energy 
T* ~ E c j « E c g 3 / 2 e-a/ 2 . The linear I-V characteris- 
tics breaks down once the voltage drop per grain exceeds 
some critical value; even in the case of the largest energy 
gap (q = 0), this value is fairly low, V c « "fE c <C E c . 

Summarizing, we have considered a Id array of metal- 
lic grains connected by highly conducting junctions. We 
have shown that the inelastic tunneling and weak charge 
quantization lead to insulating behavior below tempera- 
tures T* ~ E c e~ g / 4 (array with no background charges) 
or T* ~ E c e~ g l 2 (random array). Both scales are much 
larger than the energy E c ~ E c e~ 9 , where perturbative 
mechanisms inhibiting charge transport become sizeable. 
In essence this phenomenon is explained by the analogy 
between the array and an elastic "chain" pinned by a 
periodic potential. Most importantly, even an exponen- 
tially weak pinning potential leads to the formation of a 
"commensurate" phase where the thermodynamic DOS 
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